Scaling of the von Neumann entropy across a finite temperature phase transition 
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The spectrum of the reduced density matrix and the temperature dependence of the von Neumann 
entropy (VNE) are analytically obtained for a system of hard core bosons on a complete graph which 
exhibits a phase transition to a Bose-Einstein condensate ai T — Tc- It is demonstrated that the 
VNE undergoes a crossover from purely logarithmic at T = to purely linear in block size n 
behaviour for T > Tc- For intermediate temperatures, VNE is a sum of two contributions which are 
identified as the classical (Gibbs) and the quantum (due to entanglement) parts of the von Neumann 
entropy. 
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Introduction. Recent experimental progresses in quan- 
tum communication have increased interest on entangle- 
ment properties of quantum systems both as resource for 
quantum computingand as intrinsic characterization of 
quantum states [ll, y. A key concept to quantify the 
entanglement of a part of a system (e.g. a system of 
qubits) with respect to the rest (environment) is the von 
Neumann entropy (VNE) which has been shown to be 
related to the maximal compression rate of a quantum 
information in an ideal coding scheme for two subsys- 
tems in a joint pure state y. In the case of mixed states 
the VNE contains contributions which come both from 
classical and quantum (entanglement) correlations. In 
order to extract the entanglement part from the VNE 
the classical contribution must be subtracted, this requir- 
ing additional minimization over all possible disentangled 
states (see, e.g., [3|, |j|). The computational time needed 
to perform such a minimization grows exponentially with 
the size of the system, making the calculations impossible 
for any system with more than a few particles. 

In exceptional cases, however, the von Neumann en- 
tropy can be computed exactly due to symmetries of 
the underlying quantum problem. An example of this is 
given by integrable critical spin chains which have ground 
states exactly known in the thermodynamic limit, with 
long-distance correlations governed by conformal field 
theory [5| . In this case VNE has been proved [y, LZf to be- 
have universally as 5'(„) = | log(^ sinh(7rnT)) -|- const, 
in the limit of a large block size n and as a function of 
the temperature T, where c is the central charge of the 
underlying conformal field theory. 

In the present letter, we calculate, both analytically 
and numerically, the dependence on the temperature of 
the von Neumann entropy for a class of quantum models 
which are invariant under the action of the symmetric 
group and which arise in several physical contexts, see 
e.g. [sj, y, |lO|, |ll|, [13, [isi [iJI- A common physical prop- 
erty of these models is that they exhibit a phase tran- 
sition from an ordered to a disordered state at a finite 
critical temperature Tc [lO|. The knowledge of the ex- 



act analytical expression of the VNE together with the 
requirement that the classical part of the entropy must 
behave extensively for large n, allow us to single out 
the entanglement part of the VNE for arbitrary n and 
T. As a result, we show that the extensive part of the 
VNE coincides with the Gibbs entropy computed directly 
from the energy spectrum, while the entaglement part 
scales as S'(„) — Scibbs = 5 log n + const below Tc and as 
S(n) — Scibbs above Tc. For T = the extensive part 
Soibbs disappears and the result in Ref. [15| is recovered. 
The model. As physical model we consider a system of 
hard-core bosons described by the Hamiltonian 



H=-]^tbtb,+tb^ 



(1) 



where b\,hi, denote creation and annihilation operators 
satisfying the following hard-core Heisenberg algebra: 
[b,,bj] = [b+,b+] = Q,[b,,b+] = (1 - 2b+b,)5^j. Due to 
the on-site Fermi-like commutation relations, double oc- 
cupancy is not allowed: the action of b^ and bi on the 
single particle Fock space being 6^|0) = |1); 6i|l) = |0); 
6j|0) = 5+|l) = 0. Note that the Hamihonian ij) is in- 
variant under the action of the permutation group (sym- 
metric group Sl) and conserves the number of particles 
N = X;f=i bfbt with N < L. The ground state of the 
system in the sector with A^ particles is given by the sym- 
metric state \^{L,N)) = (^)"^^^X;p|lllllll--000000> 
where sum is taken over all possible distributions of N 
particles among L sites. Note that symmetric states are 
intensively studied with symmetry techniques for testing 
various entanglement measures, see e.g. [lg| and refer- 
ences therein. The excited states can be constructed from 
the irreducible representations (irreps) of Sl using filled 
Young tableau (YT) of type {L~r, r} with r assuming all 
values in the interval [L/2] > r > 1 where [x] denotes the 
integer part of x (see [12: llZl)- The spectral and thermo- 
dynamical properties of the system were studied in [I4I 
where it was shown that eigenvalues of ([1]) associated to 



YTs of type {L — r, r} are given by 



Er 



r+-iN{N -l)~r{r 

1j 



1)), r = 0,l, 



with degeneracy 4 = (^) - {^^i) P 



.N 



(2) 
An interesting 



property of the model is that it exhibits a phase transi- 
tion at finite temperature to a Bose-Einstein condensate 
(BEC) of hard core bosons 10|, \l^ . This can be inferred 
directly from the free energy per site F/L = K^^^J j3 
where 



in the block but only on its size n. The corresponding 
VNE was obtained in [15| where it was shown that A^ = 
(fc)(^1fc)/(^)' where k = 0, 1, ... min(n,iV). In the limit 
of large n von Neumann entropy becomes 

\ , ,1 nlL ~ n) 

Sin) « 2 l°g2 i^TTepq) + - log2 ^ . (8) 

To generalize ([8]) to arbitrary temperatures we remark 
that the reduced density matrix for temperatures T > Tc, 

L ^ cx) is 



A, 



(l-M)log(l-M)) 



min (/3^(1 - fj.) + ^ilog fi^ 

fj.E[0,min(p,qj\ 



(3) 



up to corrections of the order o{L^^). Here and below 
we denote (3 = 1/T, ^ = r/L, p — N/L and q — 1 ~ p. 
The extremum condition for Amin leads to the equation 



(3*(p* 



1 



(1-2/i* 



■In 



1 



M 



(4) 



from which we see that for T > Tc ^ {f3*{p))^^ there 
is no solution for fi* and the minimum of Q is reached 
at the end of the interval /i = min(p, g). For T < Tc, 
the minimum of ([3]) is inside the interval and is given by 
the solution of (|4]). It was shown JlCl| that this phase 
transition is actually a Bose-Einstein condensation with 
the density of particles in the condensate given by 



Pc = iP-P*{mq-P*m 



(5) 



for T <Tc and zero otherwise. In the following we shall 
characterize the behavior of the VNE across this phase 
transition. 

2. Behavior of the VNE across a classical phase transi- 
tion. At T=0 the entanglement properties of a subsystem 
of size n with respect to the rest of the system (seen as 
environment) can be characterized by the VNE 



S(„) = -tr(p(„)log2P(„)) 



;Afclog2Afc, (6) 



where A^ are the eigenvalues of the reduced density ma- 
trix p(„\ , obtained from the density matrix p of the whole 
system as /5(„) = tr(^i^_n)p. Since we are interested in the 
behavior of the VNE across the finite temperature BEC 
phase transition described above, we introduce the ther- 
mal von Neumann entropy for a block of size n as follows 

1 ^ 
S(nm = -5]rf.e-'^^'Tr(p(„)(r)log2P(„)(r)), (7) 

r-O 

where Z denotes the partition function. Note that at 
zero temperature the density matrix of the whole system 
is a projector on the completely symmetric ground state, 
p = |vE'(L, A^))(*(L,Af)| and due to the permutational 
symmetry, S(ji-j does not depend on the choice of the sites 



1 ^ 

^(") = 7 X! ^'■^"'^^'■^(n) (^) ~ P(n) U=min(p,9) , (9) 



while for T < Tc we have p(„) = p(„)(/i*) with p.* given 
by ([¥]). Thus, to compute the temperature-dependent 
von Neumann entropy we need to know the eigenvalues of 
the reduced density matrix p(„-) for arbitrary YT states. 
The following theorem provides the result. 
3. Eigenvalues of p(^n) o,nd general properties. 

Theorem. The eigenvalues of the reduced density ma- 
trix pij^\ of the eigenstates of H with N particles belonging 
to the irreps of Sl characterized by YTs of type {L~r, r}, 
with r < min{N, [i/2]) are 



A(L, N,n, r, fc, s) 



(L~n \ k — s 
\N-k 






T.^r)r--ix (10) 



k—i 



E/ ly' yjj^ \ ^ /I yn/L — N — r\/ N—r wr \ 



J=0 



Vj+i )\j+il m=C 



with fc, s, quantum numbers assuming the values k = 
0, l,...,n, and s = 0, 1, ..., min(fc, rt — fc). The correspond- 
ing degeneracies coincide with the dimension of a YTs of 
type {n,s}, i.e. deg X{L,N,n,r,k,s) = (") - (^"J. 

The theorem follows from the block diagonalization 
of the reduced density matrix with respect to the num- 
ber of particles and to the irreps of Sn. More precisely, 
P(„) can be block diagonalized with respect to the num- 
ber of bosons k appearing in the block, this leading to 
n + 1 diagonal blocks Bk, fc = 0, 1, ...,n. Each Bk can 
be further diagonalized with respect of the irreps of Sn 
which are compatible with that value of k, this leading 
to A; + 1 blocks associated to the YT of type {n — s,s} 
with s — 0, ...,min{k,n — k). Notice that the above 
decomposition implies that the block Bk has dimension 
J2i=o ((D ~ (f-i)) — ik) ^^d the dimension of the ma- 
trix p(„) is X]fc=o(fe) ~ 2", as it should be. This also 
clarifies the meaning of the quantum numbers fc, s, and 
explains the degeneracies given in the theorem. A full 
proof of the theorem will be given elsewhere. 

Before using the theorem we shall give some general 
properties of the eigenvalues (fTO]) and consider some lim- 
iting cases from which the correctness of the result ^TU\\ 
can be inferred. 



a) One can check by direct inspection that, as conse- 
quence of the above block diagonahzation, the eigenval- 
ues in pni satisfy the following remarkable sums with 
respect to s and k 



'^\ksdeg{\ks) = 



L — n \(n\ 
N-k)\k) 



s=0 
n—s 

k—s 



ik) 



ks 



/L~-n\ /L — n \ 

\r — s ) Vr+s — n— 1/ 

\r) ~ (r-l) 



(11) 



(12) 



where hereafter we use A^s as a shorthand notation for 
X{L, N,n,r,k, s). One can easily verify that the above 
sums both lead to the correct normalization of piny. 
trpin) = 1- Indeed, by using Eq. (|ll[) . we have 



trpn 



n min{k.n—k) 

E E 

fc=0 s=0 



Xk,sdeg{Xks) = ^ 



\N-k)\k) 



k=0 



(^) 



= 1. 



A similar expression is obtained by interchanging the or- 
der of the sums and using Eg. p^ . 

b) Case r = 0. In this case in the sum over m in pU]) 
only the term m = survives and, since r = implies 
s = 0, also the sum over j has only one term j = 0. The 
eigenvalues (fTOll then reduce to 



Afco — 



\N-k) 



k /L — n \ 

Sr^/k\fn-k\ yN-k) 

1=0 



(.)m = 



(^n) 



Ck) 



(13) 



reproducing the exact result obtained for the completely 
symmetric ground states in [15|. 

Although explicit (i.e. summed) expressions for the 
eigenvalues can be derived also for some other particular 
case (to be reported elsewhere), a summed expression of 
Eq. pO|) seems to be unlikely in the general case. 

In the following we will work in the thermodynamic 
limit L — > oo for which the expression for eigenvalues 
simplifies drastically. Recalling that limi^oo(''/-^) = P, 
\miL^oo{N/ L) — p, we obtain 



lim 

L — >oo 



fc) 
(^) 






hm ^(-1)' 



/N — r\/L — N — r\/r 
\j — m)\ j—m /vn 



^Ty-J, 



where 



'7 = 



(p-_m)(9_m) 
pq 



(14) 



(15) 



(16) 



plays the role of the order parameter < 77 < 1, with 
77 = for T > Tc and r; > for T < Tc. Inserting these 
expressions into (|10p and carrying out the summation 



EU vH-'^YijU) = (-^)''(l - VY we obtain 



lim Afcc 



k — s 

:p"-^'^(l-ry)^^,f('=r)("-'=-^) 

4=0 



: p'^-'^q^l ~ TiY^Fii-k + s,k-n + s-l- 7j), (17) 
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Figure 1: (a) Von Neumann entropy versus temperature 
for n — 6, obtained from Eq.© for different system sizes: 
L — 50,200,700 (squares, triangles and circles respectively). 
Continuous curves denote the corresponding quantity in the 
thermodynamic limit (|17[) while different groups of curves 
correspond to different densities of hardcore bosons in the 
system N/L = 0.1,0.2,0.5 (from down up). (b) Quan- 
tity Q(fc,s)/Qniax = Afes deg(Afes)/Qmax computed from (IT7|) 
for fixed s, versus k/n, for n = 100,200,800 (squares, cir- 
cles and triangles, respectively). Parameters are p = 0.4, 
p = 0.2 = s/n. Curves are given by the analytic prediction 
113. 



where 2^1(0, b; c; d) is Gauss hypergeometric function. In 
Figma) we compare the von Neumann entropy in ([7]) 
with the thermodynamic limit (|17p , for different temper- 
atures and various particle densities. 

4- Classical and quantum parts of the von Neumann en- 
tropy. Various limits of the above formula are discussed 
below. 

i) Zero temperature limit T —> 0. In this case rj — 1, and 
from pT|) we have A^o = p^'^q^ ('^') . This is just the hmit 
L — > 00 of Eq. (J13p , reproducing the results obtained in 

M- 

ii) High temperatures T > Tc . In this case 77 = and the 
eigenvalues (ITTl) become s-indepcndent Afc^ = p'^~^q^. 
This leads to the extensive classical entropy (Gibbs en- 
tropy) of an ideal gas with excluded volume 

Sin){T > T,) = -n{p\ogp+{l^p)\og{l~p)). (18) 

iii) Intermediate temperatures < T < Tc. This 
case < 7y < 1 is the most interesting one since 
the VNE has both classical (due to degeneracy) and 
quantum (due to entanglement) contributions. Ana- 
lyzing the sum in p7p one can show that the variable 
Q{k,s) — limi_,oo Afcs deg(Afcs) is Gaussian-distributed 
with the mean {k/n) = q, {s/n) — p: 



t^(fc, S) ~ (^max^ 



^ 2C ^ 2D 



(19) 



with Q,nax = (1 - 2p){2TTny/p{l - p){p - p){q - /i)), 
and 



C'^ = 



1 



D- 



np{l-p) n{p-p){q~pY 

(1-2a^)' ^-i {l-2p){p~q) 



n{p- p){q- pY 



B- 



n{p - p){q - l^)' 
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Figure 2: (a). Extensive part of the von- Neumann entropy 
of a block lim„^oo 5(„)/n versus renormalized temperature 
H* /p for p = 0.2. Comparison of analytic prediction (|23p 
(curve) with the exact calculations from (|17|) for finite n — 
5, 10, 50, 100, 150, 200, 300, 400 (points approaching the curve 
from above), (b). "Quantum" part of the von Neumann 
entropy S{n) — Scibbsin) versus block-size n, for p = 0.3 and 
fi* = 0.18 (points). The continuous curve refers to the fitting 
function i logj n — 0.289. 



In Figlljb) we compare the distribution of the eigenval- 
ues obtained from direct calculations with the Gaussian 
distribution p9)) . from which we see that the agreement 
is excellent. Using the above expression of Q{k, s) we 
calculate the VNE as 



n min(fc,n-fc) rifU ^ 

A:=0 s=0 



degA 



ks 



(20) 



Substituting the sums with the integrals and using the 



normalization > Q(fc, s) = 1 we obtain 



'S'(„) « -n{plog2 fJ--{l-^i)log2{l-fJ-))+-^log2n+R{q, fi) 

(21) 
where i?(q, fJ.) is a non-universal n-independent constant. 
The first term at the right-hand side of pTjl gives the 
major contribution in the large n limit, and coincides 
with the Gibbs entropy obtained from the spectrum 
of the whole system ([2|) in the thermodynamic limit 
E{L,ii,p)/Lf» H- fp +p^, dr « ,. ..l7^(^^i... , as 



Sg 



ibbs 



= lim 

L— >oc 



(3{E - F) 
L 



(22) 



where F/L — Amjn//? and Amjn is given by ^. Inserting 
the expressions for E and F into ((22|) . we obtain 

Scibbs/L = -/i* log/1* - (1 - /i*) log(l - /!*) (23) 

where fi* (T) is defined in Eq. (|4|) . The term i log2 n in 
(|2ip coincides, up to a constant, with the entanglement 
entropy at zero temperature ([S]) and can therefore be in- 
terpreted as entanglement entropy part in the von Neu- 
mann entropy at finite temperatures, defined as 



Scntin,T) = S, 



(») 



n lim 

L — ^oo 



E-F 
TL 



1 



■ log2 n 



(24) 



Note that this contribution disappears at Tc and above. 
In Fig[2]the extensive (panel (a)) and "quantum" parts 
(panel(b)) of the VNE are depicted. 

5 Conclusions. We have calculated the spectrum of the 
reduced density matrix and von Neumann entropy of a 
block of n sites, for a quantum model with permutational 
symmetry, as function of temperature and particle den- 
sity. It is shown that eigenvalues of the reduced density 
matrix in the thermodynamic limit are parametrized by 
the single parameter defined in (J16p , which also turns out 
to be the order parameter of the problem. We defined 
the entropy of entanglement for finite temperature as the 
difference between the VNE and its extensive part (co- 
inciding with the Gibbs entropy), and demonstrated its 
disappearance above the classical phase transition. For 
all temperatures below the critical one, the entropy of 
entanglement scales as ilogn. We expect this results 
to be valid also for other quantum systems with permu- 
tational symmetry (mean field models) exhibiting finite 
temperature phase transitions. 
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